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The conducting state of the quasi-two-dimensional organic conductor, Q-(BEDT-TTF)2l3, at am- 
bient pressure is investigated with ^'^C NMR measurements, which separate the local electronic 
states at three nonequivalent molecular sites (A, B, and C). The spin susceptibility and electron 
correlation effect are revealed in a locally resolved manner. While there is no remarkable site- 
dependence around room temperature, the local spin susceptibility gradually disproportionates 
among the nonequivalent sites with decreasing temperature. The disproportionation-ratio yields 
5:4:6 for A:B:C molecules at 140 K. Distinct site- and temperature-dependences are also observed 
in the Korringa ratio. Id oc (l/TiT)iK~'^ {i = A, B, and C), which is a measure of the strength and 
the type of electron correlations. The values of Id point to sizable antiferromagnetic spin correla- 
tion. We argue the present results in terms of the theoretical prediction of the peculiar site-specific 
reciprocal-space (fc-space) anisotropy on the tilted Dirac cone, and discuss the fc-dependent profiles 
of the spin susceptibility and electron correlation on the cone. 

PACS numbers: 76.60.-k, 71.20.-b, 71.30.+h 



I. INTRODUCTION 

The organic charge-transfer salt a-(BEDT-TTF)2l3 
(abbreviated as a-l^ hereafter) is a quasi-two- 
dimensional (Q2D) electron system with a 1/4-filled 
hole band which exhibits strong electron correlation 
effects [BEDT-TTF (ET hereafter) is the abbreviation of 
bisethylenedithio-tetrathiafulvalene; see Fig. 1(b)]. a-Ia 
is comprised of alternately stacked 2D conducting layers 
of (ET)+i/2 molecules and non-magnetic insulating lay- 
ers of triiodide anions (Is)^^ [Fig. 1(a)]. ^ The unit cell 
contains four ET molecules, three of which are crystal- 
lographically nonequivalent [i.e., A, B, and C molecules; 
see Fig. 1(c)]. At room temperature, the band structure 
calculation predicts semi-metallic Fermi surfaces.^ With 
decreasing temperature, resistivity shows weak temper- 
ature dependence down to Tco ^ 135 K,"^ at which the 
system undergoes a first-order phase transition from 
a paramagnetic conductor to a nonmagnetic insulator 
accompanied by an inversion symmetry breaking. ^^^"^ 
Below Tco, electrons are localized on ET lattice sites, 
forming the horizontal strip type of charge ordering 
(CO),^dO-i4 -y^fjjjch is believed to be stabilized by the 
long-range Coulomb interactions. ^^'^^ 

Under hydrostatic pressures, the phase transition 
tends to be suppressed. ^^'^''^^ Above Pq ^ 1.5 GPa, the 
charge-ordered state completely disappears and the tem- 
perature independent profile of resistivity extends down 
to several Kelvins. In this pressurized phase, recent 
band-structure calculations^^^^^ revealed that an unusual 
electronic structure with the linear dispersion is realized 
around the Fermi level, ep, (i-e., massless Dirac cone) 
which is partly confirmed by recent magnetotransport 



experiments. ^^'^'^ In contrast to the graphene monolayer, 
known as the typical massless Dirac-cone system with 
isotropic dispersion,^'* the Dirac cone in a-l^ is strongly 
tilted;^° the conical slant varies 10 times around the apex 
in the 2D fc-space.^* The tilt stems from the low lat- 
tice symmetry in a-l^ possessing only inversion symme- 
try (space group P-1). Because of this anisotropy of the 
band dispersion, a large imbalance in the local electronic 
density of states around ep is predicted among inde- 
pendent A, B, and C molecules. Recent *'^C nuclear 
magnetic resonance (^'^C NMR) experiments under hy- 
drostatic pressures observed large difference in the local 
spin susceptibilities at these molecular sites. ^^'^^ 

The Dirac-cone picture has been discussed so far only 
at high hydrostatic pressures (P > Pc) in a-l^. How- 
ever, several experiments have recently suggested that 
it might be also applicable to the low-pressure conduct- 
ing phase. Resistivity, "'27 Raman,io NMR,25'26 ^_ 
ray,^^ Ilall,^^'^^ and thermopower'^'^* measurements in- 
dicate that the ambient- and high-pressure conducting 
phases have many properties in common. Moreover, the 
band structure calculation^^ suggests that the conical 
apex (Dirac point) locates at ~ 8 meV below ep even at 
ambient pressure, ^^'^^ and can be pushed up to sp by in- 
troducing a site-dependent potential imbalance between 
A, B, and C nonequivalent molecules. ^^'^^''^^ Aside from 
this, the ambient-pressure conducting phase exhibits sev- 
eral anomalous properties such as non-Drude optical con- 
ductivity just above Tco,^^ and a large charge density 
disproportionation among A, B, and C molecular sites 
as observed in x-ray,^ Raman, ^'^ and ^^C NMR*^ exper- 
iments. On the whole, these findings probably indicate 
that the conducting state in a-l^ is closely associated 
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with the Dirac cone and/or CO fluetuations. However, 
the local electronic properties are not yet fully under- 
stood quantitatively, in particular, the electron correla- 
tion effects. 

The purpose of the present study is to clarify the elec- 
tronic nature in the ambient-pressure conducting state of 
a-Iz in microscopic and site-specific manners. Utilizing 
^^C NMR experiments, we explored the local spin suscep- 
tibility and electron correlation effects at each nonequiv- 
alcnt molecule (i = A, B, and C), since NMR is capable 
of probing the local electronic states at these indepen- 
dent molecules separately. As the hyperfine coupling of 
a nuclear spin with its surroimding electrons is highly 
anisotropic at the ^'^C site, one should know the field- 
angulear-dependent characteristics of NMR spcictra to 
make a reliable assignment of the NMR line and deter- 
mine hyperfine tensors, which correspondingly provides 
a quantitative basis for deducing the local spin suscep- 
tibility and the electron correlation effects. Thus, we 
measured the angular dependences of the NMR Knight 
shift, Ki, with rotating the external field in the planes 
parallel and perpendicular to the conducting layers [i.e., 
the ab plane in Fig. 1(c)] at each temperature measured. 
The ^'^C hypcrfinc-shift tensors and the local spin s\is- 
ceptibility were determined from the angular-dependent 
profile of Ki at each independent molecule. The nuclear 
spin-lattice relaxation rate, (1/Ti)j, was measured un- 
der the in-plane fields in order to evaluate the Korringa 
ratio, ICi oc {l/TiT)iK^'^ , which represents the strength 
and type of the electron correlations. The results are dis- 
cussed in connection with the effect of CO fluctuations 
and Dirac cone. 



II. METHODS 

To evaluate the local electronic properties in the 
ambient-pressure paramagnetic conducting state of a-Is, 
we performed ^^C NMR measurements at external field 
H of 6.00 T applied in the ah and he planes in the tem- 
perature range from room temperature down to 140 K 
(> Too)- We grew the single crystal of a-Is of the size 
2.5 X 0.5 X 0.1 mm^ by the conventional electrochemical 
oxidization method. The central double-bonded carbon 
atoms in ET were selectively enriched with ^'^C isotopes 
(nuclear spin J = 1/2, 7„/27r = 10.705 MHz/T) at 99% 
concentration [see Fig. 1(b)], which means that the ob- 
served ^^C NMR spectra come from the central carbon 
sites. The ^^C NMR spectra were obtained by the fast 
Fourier transformation of the spin echo signals. As the 
origin of the NMR shift, we referred to the resonance 
frequency of the ^^C NMR signal of TMS [tetramethyl- 
silane, (0113)481]. The spin-lattice relaxation rate, 1/Ti, 
was obtained by the standard saturation and recovery 
method. 



III. RESULTS 

A. Angular dependence of the '^^C NMR line 
shapes and line assignment 

The assignment of the observed ^^0 lines in the NMR 
spectra to A, B, and C nonequivalent molecular sites is 
the starting point of the present work. This can be ac- 
complished by considering the crystal structure, the sym- 
metry of A, B, and C molecules in the crystal, and the 
nuclear dipolar splitting in the analysis of the angular 
dependence of ^^0 NMR spectra. 

In ET-based organic conductors, all crystallographi- 
cally nonequivalent molecules contribute to distinct ^^0 
lines in the ^■^C NMR line shapes. In the conducting state 
of a-Ia, at least three peaks are thus expected associated 
with the nonequivalent A, B, and molecules [Fig. 1(c)]. 
However, each molecule has two neighboring -'^^C nuclei 
[Fig. 1(b)], which are coupled by the nuclear dipole inter- 
action. This coupling makes a single line split into a dou- 
blet or a quartet depending on whether the molecule has 
inversion symmetry or not, respectively. "^^'"^'^ The Pake 
doublet structure should be observed when the neighbor- 
ing nuclei are equivalent, which is the case in molecules B 
and 0. The Pake doublet is characterized by the nuclear 
dipolar splitting width (kHz), d, given by 

rf= 37^(1 -Ssin^C.cos^r;,), (1) 

where 7n is the gyromagnetic ratio of ^^C nucleus, 27r/i 
is Planck's constant, r = 1.360 A is the distance be- 
tween two "'^■^C nuclei,^'* and (C^, rji) are the angles formed 
between the direction of H and the ^^C=^^G vector of 
the i-site molecule {i = A, B, and C), as illustrated in 
Fig. 3(b). The quartet structure is expected when the 
neighboring ^^C sites are nonequivalent as in molecules 
A. The quartet is composed of four lines {ly = 1 — 4),^^'^^ 
the line shift (5^ of which is dependent on the nuclear 
dipolar splitting width d and the NMR shift difference 
between the two neighboring ^"^C sites, AS; i.e., Sv = 
5 ± d/3H ± v'(c^/3i?)2 + (A(5)2/2, where 5 is the aver- 
aged shift of the quartet, and equals to the average of the 
NMR shift at the neighboring ^^C sites in ET without the 
nuclear dipole coupling. Thus, in the present case, eight 
^^0 lines composed of two Pake doublets and one quartet 
are expected in the "'^^C NMR spectra. 

In order to assign the NMR lines properly, we measured 
the magnetic- field angular dependences of the NMR spec- 
tra. Figure 2 presents the typical data at 260 K with the 
external field H applied within the ab plane {H/ /ab) and 
he plane {H//be). The line shapes are strongly depen- 
dent on the direction of H. In Fig. 2(a), the ^^C lines 
can be divided into two groups, exhibiting an anti-phase 
angular dependence as the direction oi H, ijj, is changed 
within the ab plane [ip = 0° is set as H//a; see the in- 
set of Fig. 2(a)]. In Fig. 2(b), on the other hand, all the 
"'^^C lines exhibit an in-phase variation when H is rotated 
in the be plane [6 = 0° is set as H//c; see the inset of 
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Fig. 2(b)]. These phase relations can be adequately un- 
derstood from the crystal structure;"^^ i.e., the direction 
of ^'^C-2p2 orbitals within the unit cell. It is well known 
that the 2pz electrons of ^^C atoms produce a dipole- 
type hyperfine field at the ^^C nuclear position, and give 
the largest contribution to the ^^C NMR shift. ^^'^''^ As 
schematically shown in Fig. 1(c), molecules B and C, and 
molecules A form respectively a stack along the a axis. 
These stacks build a fishbone-likc molecular arrangement 
in the ab plane, and are approximately connected with 
mirror operation [{a, b, c) — >■ (a, —b, c)]. The orbitals, 
pointing perpendicular to the molecular plane [along the 
z axis in Fig. 3(a)], are parallel within the same stack, 
while they are mutually perpendicular between different 
stacks. This geometry gives rise to an anti-phase relation 
in the angular dependence of the NMR shift when the ex- 
ternal field is rotated in the ab plane. On the other hand, 
all the 2pz orbitals are roughly confined to the ab plane, 
which is expected to cause an in-phase dependence of the 
NMR lines when the filed is rotated within the be plane. 

Referring to the considerations given above, we first as- 
signed the ^^C lines in Fig. 2(a). The four lines denoted 
by open circles are attributable to the quartet (molecules 
A), and the rest four lines to two pairs of the Pake dou- 
blets (molecules B and C), as shown in Fig. 2(a). It 
is easy to distinguish molecules B and C when one fo- 
cuses on the dihedral angles between molecules A and B 
(~ 138°), and A and C (~ 131°).34 Since the former is 
larger than the latter, molecules B should exhibit larger 
phase difference from molecules A than molecules C will 
in the angular dependence of the line shapes. We thus 
assigned the ^^C lines denoted by triangles and crosses 
to molecules B and C, respectively, as shown in Fig. 2(a). 
The results arc consistent with the previous work.-'^'^ In 
Fig. 2(c), the observed angular dependences of the nu- 
clear dipolar splitting width d for A, B, and C molecules 
are presented with the corresponding symbols. The solid, 
dotted, and dashed lines are calculated angular depen- 
dences of d for A, B, and C molecules, respectively, with 
r = 1.360 A and the molecular orientations obtained from 
Ref. 34. The data agree well with the calculations. 

Next, we focus on Fig. 2(b). By comparing the NMR 
spectra at 6* = V = 90° in Figs. 2(a) and 2(b) (i.e., 
H/ /6), the ^^C lines under H/ /be were assigned as shown 
in Fig. 2(b) [with the same symbols as in Fig. 2(a)]. Fig- 
ure 2(d) presents the angular dependences of the dipolar 
splitting width d under H/ /be. Calculated angular de- 
pendences are presented together, which properly repro- 
duce the experimental results. All of these analyses thus 
assure our assignments comprehensively. 



B. Angular dependence of the ^^C NMR shift and 
determination of the hyperflne-shift tensors 

The ^^C NMR shift is proportional to the hyperfine 

field at the '^'^C site projected onto the direction of the 
external field H. The NMR shift at molecule i (= A, 



B, and C), 5^, is defined as the averaged shift of the lines 
at molecule i. In Fig. 4, we present the representative 
angular dependences of the ^'^C NMR shift 5i at 260 K 
under H//ab [Fig. 4(a)] and H//bc [Fig. 4(b)]. 6, is ex- 
pressed as 5i = [<5f^(iJ^)2 + df{Hyf + (5f (ff^)2]/|iJ|2, 
where x, y, and z arc the ET principal axes [Fig. 3(a)], 
H ={H^,Hy,H'') is the external field in the ET princi- 
pal axis [Fig. 3(b)], and Jf^, Jf^, and are the principal 
components of the ^^C-hyperfinc-shift tensor. [Note that 
the shift tensor for the molecule A with two nonequiva- 
lent ^^C sites is defined as their average, although it can 
be determined site specifically. The reason why we do so 
is that the charge and spin densities discussed below in 
terms of the shift is molecular-specific, and we need the 
shift value representative of the molecular site. Moreover, 
this definition is beneficial for the evaluation of electron 
correlation effects described in Section HID).] We fit the 
expression to the data in Figs. 4(a) and 4(b), using the 
crystal structure data given by Sawa et al.^^ The param- 
eters optimized arc listed in Table I. The fits in Fig. 4 
properly capture the sinusoidal angular dependences of 
5i under both field geometries. As seen in Table I, the 
shift tensors have different principal values among A, B, 
and C sites. This evidences differentiations in local spin 
susceptibility and charge density among these nonequiv- 
alent sites in the unit cell. 



C. Temperature dependence of the ^'^C NMR shift 
and the local spin densities 

Next, we proceed to separate the local charge den- 
sity and spin density contributions in the observed shift 
data, and investigate the disproportionation of the lo- 
cal susceptibility within the unit cell and its anoma- 
lous temperature dependence. Figures 5(a) and 5(b) 
show the temperature dependences of the ^^C NMR 
spectra at V = 110° and 50° (under H//ab), respec- 
tively. Large temperature dependences are seen in the 
NMR shift 5i when \8i\ is large. On the other hand, 
5i is almost temperature independent when \5i\ becomes 
small. The shift 5i is the sum of the Knight shift and 
the chemical shift {5i = Ki + (Ti), where the Knight 
shift, Ki, is related to the local spin susceptibility Xi 
through Ki = aiXi and the chemical shift, at, origi- 
nates from the orbital motion of electrons within ET 
molecules. '^^ is the hyperfine coupling constant be- 
tween the ^"^C nuclear and itinerant carrier spins. It is 
well known that the ^^C chemical shift in ET depends 
on the molecular valence; namely, the amount of hole on 
ET.^*" Because large charge density disproportionations 
have been observed in a-Is among nonequivalent A, B, 
and C molecules above Tco,^'^°'"'^^'^*^'^^ (Ti should be de- 
termined for each molecule in this compound. Using the 
principal values of the chemical-shift tensors af^, , 
and af^ (in ppm) determined at 60 K,^^ and linearly 
interpolating the molecular valence, pi, determined by 
X-ray diffraction^ to estimate pi at 60 K, we get the fol- 
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lowing relations between the chemical-shift tensors and 
valence p,; af^ = VSO.Spi + 46.2, af^ = 30.4p, + 159.0, 
and (jf^ = — 12.1pi + 61.9. In our case, the temperature 
dependence of ai (in ppm) is negligibly small, and the 
value of ai is determined as cta ~ 165, ctb ~ 120, and 
(Tc ~ 93 at V = 110°, and a a ~ 69, ctb ~ 165, and 
CTc ^ 166 at ^ = 50°, respectively. In general, these 
chemical-shift tensors allow us to calculate the chemical 
shift CTi under arbitrary field orientation. Thus, to de- 
termine the angular dependence of the Knight shift in 
the following, we calculated the angular dependences of 
CTi site-selectively and subtracted them from those of the 
shift 5i. 

Figures 5(c) and 5(d) depict the temperature depen- 
dences of the NMR Knight shift, Ki(= (5^ — CTj), under ip = 
110° and 50° {H/ /ah), respectively, with the use of the 
chemical shifts ct, determined above. The Knight shift 
Ki exhibits strong anisotropy within the ab plane; for in- 
stance, Ka is almost zero in Fig. 5(c), while it is large 
in Fig. 5(d). S ince the total spin susceptibility, Xspin(c'c 
2xa + Xb + Xc); is isotropic in a-Ia,^ the anisotropic be- 
havior of Ki should stem from the anisotropy of the hy- 
perfine interaction at the ^^C position, which is expressed 
as Ki = [af'{H'=)^ + af (T??')^ + al''{H^f]xi/\H\^ with 
af^, , and af^ the principal components of the ^^C 
hyperfine-coupling tensor [sec Fig. 3(a)]. [This means 
that the small value of K^ in Fig. 5(c) should be at- 
tributed to the vanishingly small hyperfine-coupling con- 
stant at ij) = 110° geometry.] In ET-based materials, it is 
known that the profile of the highest occupied molecular 
orbital (HOMO) of ET determines the anisotropy of the 
coupling tensor af^ {ji = x, y, and z).^^ The temper- 
ature dependence of is negligible, since the spatial 
distribution of HOMO is most likely to be temperature 
independent. Moreover, the spatial profile of HOMO 
is reasonably assumed to be the same at all molecules 
in the unit cell because the observed differences in the 
molecular structures are very small. '^^ Therefore, in the 
first approximation, we can assume that a-*'* is site- and 
temperature-independent. We also note that molecules 
A, B, and C arc arranged in a nearly symmetrical manner 
within the ab plane, as we mentioned above [Fig. 1(c)]. 
Hence, if H is rotated in the ab plane, the average, i^f™, 
and the amplitude, K^™"^, of the angular dependences of 
the Knight shift [see Fig. 6] are both expected to reflect 
the magnitude of xf, ^f™ = ^^''X^ and K-""^ = a^'^Pxi 
with a**™ and a**™? the averaged hyperfine-coupling con- 
stants. 

Figure 7 shows the temperature dependences of Kf^'^ 
and K^™^. For determining them, we calculated the an- 
gular dependence of ai, using the empirical relations de- 
duced above, and subtracted them from the angular de- 
pendences of the shift, 6i, at all temperatures measured. 
The external field H was rotated in the ab plane, because 
iff™'' becomes largest in this field geometry and the line 
shapes are readily discerned owing to the anti-phase an- 
gular dependence of the lines as mentioned above (in 
Sec. Ill A). Around room temperature, there are lit- 



tle site-dependences in the observed iff™ and iff™''. 
With temperature decreased, however, they exhibit site- 
specific temperature dependences, and then begin to de- 
crease monotonically below T ~ 180 K at all sites down 
to Too ~ 135 K. 

To deduce the local spin susceptibility Xi from iff™ 
(iff'"P)^ it is necessary to evaluate the hyperfine-coupling 
constant a^™ (a'^^P) defined above. This is achieved 
by comparing iff™ (iff™'') with bulk spin susceptibility 
Xspin at room temperature where the site dependence of 
^avo (xf"iP) becomes negligible; that is, iff™ oc Xspin 
(iff™'' cx Xspin)- At room temperature, iff™ w 360 
and i^f™P « 440 (in ppm) with little site dependences 
(Fig. 7) and Xspin = 6.8 x 10~* emu/mol f.u. mea- 
sured by Rothacmcl et al.^ In terms of them, the hy- 
perfine coupling constants are evaluated as a'^'^'^ k, 5.9 
kOe//iB and a"*™ w 7.2 kOe//XB. To check whether 
these coupling constants are applicable under arbitral 
temperatures, we compare the temperature dependences 
of iff™/4 (sharps) and Y.i Kf^^ (squares) to that 
of the spin susceptibility Xspin (crosses), as depicted in 
Fig. 8, where the summation is taken over all molecules 
in the unit cell (i.e., molecules A, A, B, and C; see Fig. 1). 
All of these quantities are nearly scaled to each other (i.e., 
Ei^f"74 <x Xspin and Y.iK""^!'^ ^ Xspin), which as- 
sure our assumptions that a^™ and a^^'^ are site- and 
temperature-insensitive. Thus, it is allowed to determine 
the local susceptibility x% with these coupling constants 
through Kf"" = a^™Xi and iff™'' = a'*"Pxi at all tem- 
peratures measured. 

In Fig. 9(a), we show the temperature dependences 
of the hereby deduced local spin susceptibility Xi from 
Figs. 7(a) and 7(b). [Note that Xi in Fig. 9(a) is an 
average among iff™/a'*™ and /ff™P/aamp ijgcause these 
values are in good agreements with each other over the 
entire temperatures measured.] Susceptibility shows lit- 
tle site dependence aroimd room temperature. With de- 
creasing temperature, however, large imbalance develops 
in Xi among A, B, and C nonequivalent molecules below 
^ 270 K, which increases down to Tco- Below ^ 180 
K, all Xi's decrease monotonically with decreasing tem- 
perature. The tendency, xc > Xa > Xb, is consistent 
with the preceding works. ^'^^^^'^^ Just above Tco, its ra- 
tio reaches 5:4:6 for A:B:C molecules. This large imbal- 
ance cannot be expected by a simple single-band picture. 
Notice that the temperature dependence of XA has the 
same profile as that of Xspinj namely, XA oc Xspin, as seen 
in Fig. 8 (circles). 

To highlight the temperature dependence of the spin- 
density disproportionation prominently, we depict the 
relative local spin density in Fig. 9(b), defined as (sj) = 
Xi/[(2xA + Xb + Xc)/4], which reflects the relative con- 
tribution of the i-site HOMO to the conduction band 
around ep- Imbalance develops in the relative local spin 
densities at the B and C sites with decreasing tempera- 
ture. On the other hand, the A site shows little temper- 
ature dependence ((sa) ~ 1), which is the direct conse- 
quence of Xspin C!C XA- At first glance, it seems to be plau- 
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siblc that the observed spin-density disproportionation 
can bo quahtatively understood with the semi-metalhc 
band picture as a redistribution process of B- and C-site 
HOMO contributions in the two bands. Indeed, band- 
structure calculations predict two bands locating in the 
vicinity of the Fermi level £f-^^ However, it is difficult to 
explain the strong decrease in Xspin (Fig- 8) by the sim- 
ple semi-metallic framework; e.g., in the case of 2D Fermi 
surfaces, Xspin[oc D{eF)] is expected to show little tem- 
perature dependence, since the density of states D{ep) is 
not sensitive to the value of sp. Moreover, the trend of 
the spin-density disproportionation is opposite from what 
is expected from the x-ray diffraction measurement^ [B 
site: (hole) rich, and C site: (hole) poor], which is also 
anomalous as a conventional semi-metal. As we shall see 
in Sec. IV A, these features should be rather regarded 
as consequences of the anisotropic Dirac cone realized 
around ep, which was originally predicted under hydro- 
static pressures in a- 13.^° 

D. Temperature dependence of the ^'^C nucleEir 
spin-lattice relaxation rate and the electron 
correlation effects 

So far, the local electronic states in a-Is have been 
revealed in terms of the static susceptibilities. The nu- 
clear spin-lattice relaxation rate, l/Ti, probes the fluctu- 
ations of electron spins. In this section, we show the site- 
dependent spin dynamics uncovered by the site-selective 
measurements of 1/T\T. which allows one to see the cor- 
relation effects in site-specific or band-specific manners, 
as discussed in Sec. IV B. 

In the B and C molecules with inversion center, the two 
neighboring carbons are equivalent and give the identical 
relaxation rate. However, the two nonequivalent carbons 
in the molecule A without inversion center should exhibit 
different relaxation rates due to different hyperfine cou- 
pling constants as reported earlier.''^ Nevertheless, the 
so-called T2 process, which works to average the site- 
specific relaxation rates among these two carbons, tends 
to alter the two values toward some intermediate val- 
ues in-between them. Then, the distinction of the two 
observed rates is not so informative, but their average 



is a meaningful value specific to the A molecule. Thus, 
we determined the relaxation rate representative of the 
molecule A from the relaxation curves of the whole spec- 
tra (i.e., the quartet), which were nearly single exponen- 
tial. 

In Fig. 10(a), we present the temperature dependence 
of the ^^C nuclear spin-lattice relaxation rate, 1/Ti, di- 
vided by temperature T at 7/; = 110° imdcr H/ /ah. With 
decreasing temperature, 1/TiT decreases monotonically 
at all sites from room temperature down to Tqq. The 
temperature dependence is, however, different from site 
to site. The decreasing rates are large at A and B sites, 
while C site exhibits only moderate temperature depen- 
dence. In a conventional single-band metal, 1/TiT is not 
expected to show site-specific temperature dependences 
because all sites probe the same electronic properties in 
the conduction band. The distinct behaviors at A, B, 
and C sites indicate that the electronic structure in this 
system cannot be interpreted within a simple single-band 
model. 

Since the values of l/TiTK"^ measures the degree of 
electron correlations in the conducting state, l/TiTK"^ 
was evaluated at each site. However, the Knight shift 
at A site is too small ( if a = —60 ~ —70 ppm) at this 
filed geometry in Fig. 5(c) (■)/' = 110°) to obtain reli- 
able values of 1/TiTK^, compared to B and C sites. Wc 
thus measured the spin-lattice relaxation rate 1/Ti on 
A site at V = 50°, where the Knight shift at A site is 
large and hence the relative error in 1/TiTK^ becomes 
small. The results are shown in Fig. 10(b). Then, we 
evaluated the values of l/TiTK"^ for the data series at 
ip = 50° for A site, and at ijj = 110° for B and C sites, 
respectively. In the present case of anisotropic hyper- 
fine couplings, and in the presence of electron correla- 
tions, 1/TiTK'^ is expressed as the modified Korringa 
relation; {l/TiT\K-^ = (47rfcB/ft)(7n/7c)'ft(C*, ^^O-^i (« 
= A, B, and C).''^''''^ Here. 70 is the gyromagnetic ratio of 
an electron, is the Boltzmann constant [which yield 
(47rA:B/?i)(7n/7e)' = 2.397 x lO^ sec"! K"!], /3.(C.,r;.) 
is the i-site correction factor for the anisotropy of the 
hyperfine-coupling tensor as defined in the following (i = 
A, B, and C) 



(af^/af^)^(sin^ rji + cos^ Ci cos^ rji) + (a|'^/af^)^(sin^ r]i + cos^ Ci cos^ r]i) + sin^ ^ 
2[«^/af^) sin^ Ci cos2 rji + (af /a?^) sin^ Ci sin^ m + cos2 Ci]^ 



(2) 



and /Ci is the i-site Korringa ratio, which reflects the type 
and strength of the electron (;orrdations. The direction 
of H, {Ci,rii) [see Fig. 3(b)], is determined at i = A, 
B, and C molecules for the present field geometries, as 
shown in Table II, based on the molecular orientations in 



Ref. 34. The principal components of the ^''C hyperfine 
coupling tensor, af^, af^, and af^, are determined at each 
molecule i (listed in Table II) from the total shift tensors, 
given in Table I, and the estimated chemical-shift tensors 
at 260 K (see Sec. IIIC). The values of /3i(Ci,?7i) are also 
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presented in Table II. 

Figure 11 represents the temperature dependence of 
the Korringa ratio /Ci[oc {l/TiT)iK^'^]. JCi = 1 means 
that there is no electron correlations. Figure 11 shows 
relatively large JCi for all sites (/Cj ~ 4 — 11), indicat- 
ing the presence of antiferromagnetic correlations. At 
room temperature, ICi shows little site dependence. Upon 
cooling, however, it begins to exhibit distinct behaviors 
for A, B, and C molecules below ^ 270 K. The Kor- 
ringa ratios exhibit only small temperature dependences 
at A and C sites (/Ca ~ 7 and /Cc « 4 — 6), while 
/Cb increases below 200 K and reaches a value of 11 
at 140 K. In the lowest temperature region, the mag- 
nitude relation of /C, is given by Kb > Ka > )Cc ■ 
We note that the well-studied dimer-type organic metal 
K-(BEDT-TTF)2Cu[N(CN)2]Br,33.35'40 ^j^j^jj resides on 

the verge of the Mott transition, shows /C ~ 5 — 10. On 
the other hand, in 6l-(BEDT-TTF)2l3,''^ known as a good 
2D metal with a 1/4-fillcd band, exhibits K, « 2—3, point- 
ing to weak correlations. In our case, JCi is intermediate 
or as large as in K-type salt. This indicates the impor- 
tance of electron correlations in a-Ia, which is discussed 
in more details in Section IV B. 



IV. DISCUSSION 

As we mentioned in the preceding sections (Sec- 
tions IIIC and HID), it is difficult to explain the ob- 
served anomalous decreases in spin susceptibility Xspin 
(Fig. 7), temperature-dependent spin density dispropor- 
tionations (Fig. 9), and the large temperature and site 
dependences of spin fluctuations /Cj[oc {l/TiT)iK^'^] 
(Fig. 11) in terms of simple semi- metallic or single band 
pictures. Then, how can we understand our NMR re- 
sults comprehensively? Intriguingly, there are two hints 
in the previous works. First, as we noted in Sec. I, re- 
sistivity shows weak temperature dependence in this sys- 
tem both at ambient and high pressures. Secondly, 
the NMR local spin susceptibilities at intermediate^^ and 
higli-^' pressures exhibit very close behaviors to our re- 
sults at ambient pressure. These results suggest that 
the ambient- and high-pressure conducting states have 
qualitatively similar features. Furthermore, the observed 
spin-density disproportionation (Fig. 9) reflects the pres- 
ence of site-dependent potentials. The trend of the dis- 
proportionation is the same range as that expected in 
the local site-potentials in Ref. 18, which is predicted 
to stabilize the Dirac cone.^^'^°'^^'^° All of these consid- 
erations thus imply that the Dirac cone dominates the 
electronic properties even at ambient pressure. In fact, 
as shall be discussed below, the observed features are 
properly captured by the theoretical consequences of the 
anisotropic conical dispersion realized around ep-"^^ We 
will focus on the local spin density disproportionations in 
Sec. IV A, and the site-dependent spin fluctuation eff'ects 
in Sec. IV B. 



A. Disproportionation of the local spin densities 
and the site-specific fc-space anisotropy in the 
conduction band 

First, we consider the temperature and site depen- 
dences of the local spin densities. Let's assume that the 
wave function in the conduction band is given by the 
Bloch sum of the HOMO (highest occupied molecular 
orbital) at the z-sitc ET molecule {i = A, B, and C), 
(j)i[= (j){r — Ri — Si)], as follows: 

*feW=E E e"'-'''Ci,k<f>ir-Ri-Si), (3) 

Rl i=A,B,C 

where Ri stands for the position of one a-l^ unit cell, and 
Si is the vector connecting the z-site to the A site in the 
unit cell {Sa = 0). The experimentally obtained i-site 
local spin susceptibility Xi is proportional to the thermal 
average of 4'i contributing to the conduction band around 
^F, (Efe |Ci,fcP>^=^^±fe3T, which is explicitly given as^^ 
X. = -/_°LrfeA(e)/'(e) = - rfe 1^,^!^% - 
^fe)/'(^)> with /(e) the Fermi-Dirac distribution function, 
^fc the energy-momentum dispersion of the conduction 
band, and Di{e) the local electronic density of states at 
molecule i. Prime stands for the derivative with respect 
to s. 

There are two significant consequences of the band- 
structure calculations relevant to the present results;^^ 
one is the large tilting efifect of the Dirac cone, and 
the other is the resulting strong angular dependences of 
IC'i, fell's about the cone. The ratio of the steepest and 
gentlest slants (i.e., anisotropy of the Fermi velocities) is 
estimated at about 10,^^ and the latter leads to a flat 
band dispersion giving the van Hove singularity around 
10 meV above the Dirac points, ^-'^ as schematically de- 
picted in Fig. 12. Noticeably, |CB,feP is largest around 
the steepest dispersion, denoted as S in Fig. 11, and 
shows a node around the gentle dispersion, denoted as 
G, whereas |Cc,fcP has opposite characteristics. |CA,fc|^ 
is predicted to show no remarkable angular dependence. 
This means that the local spin susceptibilities at B and C 
sites preferentially probe the thermal excitations around 
the S and G portions of the cone, respectively, while the 
susceptibility at the A site sees the average over the cone. 
The overall temperature dependence of the spin suscep- 
tibility, shown in Fig. 8, is basically explained in this 
context as follows: the conical dispersion has an energy- 
linear density of states, which gives rise to a linearly 
temperature-dependent spin susceptibility. ^^'^^ The de- 
creasing total susceptibility with lowering temperature, 
observed below 200 K in Fig. 8, can be thought of as signi- 
fying this. The leveling-off of the susceptibility at higher 
temperatures, on the other hand, implies the breakdown 
of the cone picture at high energies. Actually, the calcu- 
lated total spin susceptibility based on the band structure 
tends to saturate at high temperatures.^^ 

The spin density disproportionation of xb < Xc, en- 
hanced at low temperatures [see Figs. 9(a) and 9(b)], 
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reasonably corresponds to the small and large local den- 
sity of states at the S and G portions on the cone. The 
peak formation of xc around 190 K is attributable to 
the van Hove singularity located around G, where |Cc,fcP 
shows the maximum. On the contrary, |CB,fcP is van- 
ishingly small at G, as mentioned above, which explains 
why xb continues to increase monotonously up to room 
temperature. Meanwhile A site shows temperature de- 
pendences in between B's and C's both in the local spin 
susceptibility XA and the relative local spin density (sa) 
[see Figs. 9(a) and 9(b)]. xa scales to the bulk spin 
susceptibility, Xspin, over the whole temperature range 
(xa oc Xspin; see Fig. 7), which is reflected in the tem- 
perature independent profile of (sa) shown in Fig. 9(b). 
These results indicate that A site probes the whole con- 
ical dispersion in an averaged manner, and supports the 
theoretical prediction that |CA,fcP shows no remarkable 
angular dependence on the cone.^'^ 

At high temperatures of hundreds Kelvin, the elec- 
tronic mean free path in molecular conductors can reach 
the order of the unit cell size due to the strong electron- 
phonon scatterings. The high-temperature equalization 
of the local spin susceptibilities [sec Fig. 9(a)] might be 
partially due to this scattering effect, which will ther- 
mally average the anisotropy of the wave functions in 
the fc-space. Several organic conductors exhibit bad 
metal natures such as a loss in Drude weight in optical 
conductivity, or broadening in ARPES spectra at high 
temperatures.''^'^'* Moreover, the correlation-induced ef- 
fect such as charge ordering is disturbed at high tem- 
peratures in general. In the present system, the charge 
disproportionation, if it is enhanced by electron correla- 
tions, is expected to be depressed at high temperatures. 
This effect tends to make the Dirac point pushed down 
below ep, and may lead to a crossover from the Dirac 
system to a semi-metal. Actually, the spin susceptibil- 
ity and Korringa ratio arc both temperature- and site- 
independent around room temperature. It is probable 
that the picture of the massless Dirac electrons is broken 
by these thermal effects. 

To illustrate how the hole density at the molecule ? , pi, 
relates to the above discussed i-site spin susceptibility Xi, 
we made a comparison between pi and Xi ■ Here, the value 
of Pi is estimated from the intra-molecular bond lengths 
determined by the X-ray diffraction study by Kakiuchi 
et al.^ and the charge-sensitive modes in the vibrational 
spectroscopy by Wojciechowski et al.^'^ They found that 
Pb > PC I which is supported by Hartree-Fock calcula- 
tions based on the extended Hubbard model. This re- 
lation is, however, opposite to the spin density profile, 
Xb < Xc- Because our material is a 3/4-filled electron- 
band (or a 1/4-filled hole-band) system (see Sec. I), the 
charge (hole) density might correlate to the spin den- 
sity, which seems irreconcilable with the present results at 
first glance. However, the charge-spin correlation holds 
when both charges and spins are spatially well local- 
ized. For itinerant electron systems, on the other hand, 
they should not necessarily match since the amount 



of the hole is given by pi cx X]fc|C'«,feP[l - .f(Cfe)] = 
(Efe |Ci,feP)e>£^ [with J2iPi = 2 («: aU molecules 
in unit cell)], while the spin density is expressed as 
X» oc {J2k\'^hk\'^)^=^j,±kBT mentioned above. Ob- 
viously, Xi is determined by jCi^fej^ only in the vicinity 
of e_F, while pi reflects the whole summation of jCi^fej^ 
over the conduction band above sp. Hence, we acquire 
the following relations from the experimental results: 
(EJC^B..|^),^,,±,,^ < (EJCc..!^),^,^^,^^ and 
(EJC'B,fe|^),>,^ > <Efe|C^c,fe|'),>,^, which suggest 
that Efcl^B.fej^ < Efc|C'c,fcP holds near sf whereas 
Efe |C'B,fcP > Efe |C'c,fcP is valid weU above ep- These 
contrasting relations at low and high energies are consis- 
tent with the picture of the conical band dispersion which 
is characterized by the distinct anisotropics in |CB,fc|^ 
and |Cc,fc|^;^^ the portion G (in Fig. 12), which is con- 
tributed largely from the site C, is cut off by the Brillouin 
zone boundary at low energies, while the portion 8 rele- 
vant to the site B extends to higher energies and over a 
wide range of the fe-space covering the zone center (see 
Fig. 12). 

B. Korringa ratios and spin fluctuation effects 

Next, we turn our attentions to the site-specific spin 
fluctuation effects. As noted in Section HID, the Ko- 
rringa ratio, JCi oc {l/TiT)iK^'^ , is a measure of the 
strength and type of spin fluctuations in the conduction 
band.^^ As represented in Fig. 11, /Cj is strongly site 
and temperature dependent in our system, in contrast 
to a simple metal where ICi should be independent on 
sites and temperatures. This anomalous behavior is un- 
derstood based on the theoretical consequence that each 
molecular orbital at the site i contributes to the coni- 
cal dispersion with a differing angular dependence in the 
fc-space, which is consistent wi th the discussion in the 
preceding section. 

First, it should be reminded that |CA,fc|^ has a little 
angular dependence in the reciprocal space. Hence, ICa 
is expected to probe the fluctuations over the entire Dirac 
cone on average and to serve as a benchmark for the 
strength of correlations in the system. As we mentioned 
above, /Ca exhibits a relatively large value (/Ca ~ 7) 
in the whole temperature range (Fig. 11), pointing the 
presence of strong or intermediate antiferromagnetic spin 
fluctuations in the system. 

The spin-lattice relaxation rate at the site A, (l/ri)A, 
is given by*^ 

(l/Ti)A = 2^-j -^(aA) 2.^^> (4) 
' Q 

where is the transverse component of the hyperfine 
coupling tensor at the site A, to is the NMR resonance 

frequency, and x"{Qi'^) is the imaginary part of the dy- 
namical spin susceptibility. The wave vector, Q, and 
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the frequency lj are characteristic for the spin fluctua- 
tions associated with the electron-hole pair excitations 
a,i ep- For the free electrons, x"(Q,a;) is expressed as 
XoiQ,^) = 7r72;i3u;[Aot(£F)]/2,4^ with Aot(eF) the to- 
tal density of states at £:_f. In the presence of spin fluctua- 
tions, on the other hand, x"(Q; £^"6 enhanced from the 
constant value Xo(Q,w) at the wave vector Q connect- 
ing the degenerate states around ep- The present system 
has two Dirac cones (at valley fco and — fco) in the first 
Brillouin zone, which are mutually connected with time- 
reversal symmetry.^*^ {k-Q stands for the position of the 
Dirac point inside the first Brillouin zone.) In this situ- 
ation, two kinds of interactions arc conceivable: (i) the 
intra-valley scattering with Q ^ 0, and (ii) the inter- 
valley scattering with Q ~ 2feD- The former scattering 
gives ferromagnetic fluctuations while the latter leads to 
antiferromagnetic ones. The result of /Ca « 7 in much 
excess of unity clearly shows that the inter- valley scatter- 
ing is the dominant source of the enhanced spin fluctua- 
tions in our system. This enhanced electron correlation is 
naturally understood to help stabilizing the correlation- 
induced CO-insulating state below 135 K.^^'^^ 

In contrast to |CA.fc|^, IC'e.fel^ and |Cc.feP have so dis- 
tinct angular and energy dependences on the cone that 
/Cb and /Cc are expected to measure the fluctuation ef- 
fects in preferential areas of the fe-space; namely, the S 
and G portions in Fig. 12, correspondingly. Note that, 
roughly speaking. Id probes the intensity of spin fluctua- 
tions per thermally excited quasi-particles instead of the 
intensity itself because (l/TiT)j is divided by Kf. At 
low temperatures, /Cb is much enhanced compared with 
/Cc; at 140 K, the value of /Cb reaches 11, which is twice 
of that of /Cc- The thermally excited states probed at 
B site are well located around the Dirac point in the 2D 
fc-space because the dispersion of the S portion is steep 
(see Fig. 12). In the G portion, on the other hand, the 
thermal excitation should be extended in a wide range of 
fc-states because of the flat dispersion with the van Hove 
singularity located around 150 K above ep- The different 
fe-space proflle of thermal excitations probed at B and C 
sites should give a corresponding difference in the inter- 
valley scattering seen at B and C sites. The inter-valley 
scattering probed at B site are well restricted to the vicin- 
ity of Q ~ 2feD, whereas that probed at C site should be 
spread aroimd Q ^ 2feD- The present results imply that 
the antiferromagnetic fluctuations are sharply enhanced 
at a wave vector of Q ~ 2fcD, and that the inter- valley 
scattering is responsible for the electron correlations in 
this system. 

At elevated temperatures, /Cb and /Cc gradually ap- 
proach each other, and all /C^'s become nearly the same 
around room temperature (/C^ « 6). This temperature 
dependence is understood in line with the diminishing 
spin density disproportionation at higher temperatures, 
as discussed in Sec. IV A [see Fig. 9(b)]; the fe-space av- 
eraging due to the intense phonon excitations makes the 



difference between the molecules A, B, and C indistin- 
guishable, and/or the possible depression of charge dis- 
proportionation pushes down the Dirac point below ep, 
causing a semi-metallic state. 

Theoretically, electron correlations stabilize the Dirac 
cone dispersion in this system as we mentioned in 
Sec. 1.18,20,29,30 Furthermore, the presence of the tilted 
Dirac cone can largely enhance the anisotropy of |CB,fcP 
and |Cc,fe|^ in the fe-spacc.^-'^ The present site-dependent 
NMR characteristics well below room temperature are 
consistent with these predictions and give an opportu- 
nity to look into the magnetic properties of the tilted 
Dirac cone in a fe-dependent manner. 



V. CONCLUSION 

Wo performed ""^^C NMR measurements on the charge- 
disproportionated conducting state in the layered organic 
conductor Q;-(BEDT-TTF)2l3 at ambient pressure. Re- 
flecting the low crystal symmetry, we obtained separate 
NMR lines for crystallographically nonequivalent three 
molecules in the unit cell {i = A, B, and C). The tem- 
perature dependences of the resultant site-specific Knight 
shift Ki are properly captured by the conical disper- 
sion with the Dirac points close to £p. The analyses 
of the site-selected nuclear relaxation rate (l/TiT)j and 
Ki point to the presence of strong or intermediate anti- 
ferromagnetic spin correlations. Exploiting the theoreti- 
cal prediction of the peculiar site-specific reciprocal-space 
anisotropy in the Dirac cone, the present results turn out 
to show that the local spin susceptibility and electron 
correlations are both strongly angular dependent on the 
cone. This outcome is regarded as one of the outstanding 
aspects inherited from the tilted Dirac cone in Q!-(BEDT- 
TTF)2l3. 
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(a) I3 (b) (c) 




FIG. 1: (Color online) (a) Side view of the crystal structure 
of Qf-Is- (b) Molecular structure of BEDT-TTF (ET). Arrows 
indicate the position of carbon atoms enriched by ^^C isotopes 
with 99% concentration, (c) Top view of the schematic crys- 
tal structure in the conducting layer. Rectangular represents 
nonequivalent ET molecules {— A, B, and C) in the unit cell. 
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TABLE 1: The /i-axis component of the ^"^C-hyperfine-shift 
tensor 5f {fi = x, y, z; see Fig. 3) at molecule i { — A, B, 
and C) deduced from the angular dependences of the NMR 
shift Si's at 260 K shown in Figs. 4(a) and 4(b). 
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FIG. 2: (Color online) (a), (b) Angular dependences of the ^^C NMR line shapes under (a) H//ab and (b) H//bc at 260 
K. The definitions of angles ip and 9 are given in the insets. [Note that the origins of tp and 6 are set as H//b and H/fc, 
correspondingly], (c), (d) Nuclear dipolar splitting width d at A, B, and C nonequivalent molecules extracted from Figs. 2(a) 
and 2(b), respectively. Calculated angular dependences are shown together based on Ref. 34 (curves). In all figures, symbols 
represents the data relevant to molecules A (open circles), B (triangles) and C (crosses), respectively. 









af 








i 


(kOe//iB) 


(kOe/MB) 


(kOe//iB) 


(degree) 


(degree) 




A 


-1.11 


-1.61 


13.35 


21.4 


57.2 


0.106 


B 


-1.26 


-2.25 


11.91 


42.6 


67.7 


1.139 


C 


-0.74 


-0.70 


14.32 


31.8 


58.8 


0.281 



TABLE II: Principal components of the hyperfine-coupling tensor af^ = x, y, z) at molecule i, spherical polar angles {C,i, 
Tji) [defined in Fig. 3(b)], and the correction factors for the anisotropy of the hyperfine-coupling tensor l3i{C,i,'r]i) [see Eq. (2) in 
Sec. HID] at i = A (V- = 50°), B (?/> = 110°) and C (V' = 110°) molecular sites. 
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FIG. 3: (Color online) (a) Molecular principal axes (a;, y, z) 
of ET molecule, (b) Spherical polar coordinate (C,i, rji) of the 
external field vector H at i-site ET molecule (i = A, B, and 
C) in the principal axes {x, y, 2).^^ 
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260 K, he plane 
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6 (degree) 

FIG. 4: (Color online) Angular dependence of the ^^C NMR 
central shift at molecule i (= A, B, and C), 5i, under (a) 
H//ab and (b) H//bc at 260 K. Open circles, triangles and 
crosses stand for the nonequivalent molecules A, B, and C, 
respectively. The solid, dotted, and dashed lines are fitted 
curves 5, = [5r{H^f + S^''{Hyf + Sr{H'f]/\H\^ (see text) 
with the optimized parameters of Sf^, 5^^ , and J" listed in 
Table I. 
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FIG. 5: (Color online) (a), (b) Temperature dependences of the ^'^C NMR spectra under H applied within the ab plane {H/ /ah) 
at (a) i/; = 110° and (b) V = 50°. [The angle V is defined in the inset of Fig. 2(a).] (c), (d) Corresponding "C NMR Knight 
shifts at A, B, and C nonequivalent molecules. 




FIG. 6; (Color online) Definitions of the average and ampli- 
tude of the angular dependence of the ^'^C NMR Knight shift, 
Kr" and 
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FIG. 7: (Color online) (a) Average (J^r") and (b) amplitude 
^^amp-j ^j. ^j^g sinusoidal angular dependence of the ^^C NMR 
Knight shift at i-site ET molecule {i = A, B, and C) under 
the external field H applied parallel to the ab plane {H/ /ab). 
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FIG. 8: (Color online) Temperature dependences of the nor- 
malized susceptibility Xspin (crosses), from Ref. 8, Kf" /A 
(sharps), /4 (squares), and xa (circles), respectively. 

Summations are taken over all molecules in the unit cell. 
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FIG. 9; (Color online) (a) Temperature dependences of the 
i-site local electron spin susceptibility, Xi (* = B, and 
C nonequivalent molecules), (b) Relative local spin density, 
=X»/[(2XA + XB + Xc)/4]. 
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FIG. 10: (Color online) "C NMR nuclear spin-lattice relax- 
ation rate divided by temperature, 1/T\T, at A, B, and C 
nonequivalent molecules under (a) i/; = 110° and (b) ^jj = 50° 
{H//ab). 
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FIG. 11: (Color online) "C NMR Korringa ratio K.i[(x 
{l/TiT)iK~'^] at i = A, B, and C nonequivalent molecules. 
Dashed line {Id — 1) corresponds to the free electrons case. 
)Ci > 1 implies the antiferromagnetic correlations, while 
/Ci < 1 stands for the ferromagnetic ones. 
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FIG. 12: (Color online) Schematic representation of the fe- 
space anisotropy of the tilted Dirac cone spectrum predicted 
by the band-structure calculation given in Ref. 21. S and G 
stand for the steepest and gentlest portions of the dispersion, 
and feD ( — ^d) represents the position of the Dirac point. (In- 
set) First Brillouin zone of a-Ia. Inter- valley scattering vector 
with the wave vector Q ~ 2ko is represented by the bold ar- 



row. 



